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Abstract 

We prove spectral and dynamical localization on a cubic-lattice quantum 
graph with a random potential. We use multiscale analysis and show how 
to obtain the necessary estimates in analogy to the well-studied case of 
random Schrodinger operators. 



1 Introduction 

Since the middle of the 1980's the mathematical approach to the phenomenon of 
localization in random solids witnessed a rapid development. One of the tech- 
niques used to prove localization is multiscale analysis. Introduced by Frohlich 
and Spencer in [FS83] and further developed by von Dreifus and Klein in [DK89J 
for the original Anderson model on the lattice, it had been extended to the con- 
tinuum by Combes and Hislop in |CH 94]. By now there is a large number of 
discrete and continuum models for which localization has been established this 
way, see |Sto01| and for more recent advances [GK01J. 



On the other hand in recent years the interest also turned to the shape of struc- 
tures made of semiconductor and other materials. In particular, quantum graph 
models became popular as models of various superlattice structures. Therefore it 
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seems natural to ask how one can extend the multiscale proof of localization to 
such graph models. In this paper we want to give an answer for a particular case 
of a special cubic lattice graph that can be embedded in R d , so that the known 
techniques work similarly as in the "continuum" case. Recall that rectangular 
lattice graphs also exhibit other interesting spectral properties |Ex95| . 

The embedding into R d provides an easy way to describe our graph T. Let 
V(r) = Z d be the vertex set and let the set of edges E(T) consist of all line 
segments of length one between two neighbouring vertices in directions of the 
coordinate axes. As usual we identify each edge with the interval [0, 1] with ori- 
entation in the sense of the increasing coordinate in M. d . The initial and endpoint 
of an edge e are labeled by t(e) and r(e). 

The embedding of T into M d allows us to define subgraphs of T in terms 
of suitable domains in M. d . To make this precise, we will call a bounded domain 
A C M d with piecewise smooth boundary T—edge bounded (T-ebdd.) if <9A C E(T) 
and for each edge e G E(T) either e C <9A, or e intersects dA at most in its 
endpoints. The graph T fl A arises from T by deleting all the edges outside A 
(including those on the boundary). For its sets of edges and vertices we write 
E(T H A) and V(r n A), respectively. 

The Hilbert space underlying our model is ^(r) := © e& g(n -^2(0, 1); in a 
similar way we associate L 2 (T n A) := © ee £( rnA ) £2(0, 1) with T fl A. Further we 
need the Sobolev space of order one, 

Wl(T) := {/G0 W%(0, 1) I / continuous at all vertices v G V, 

WfWw^T) ■= Yl 11^11^(0,1) < °°} 

ee£(r) 

with the obvious notation and terminology for edge components of /, and its 
analogue W^iV n A). 

We can now define the random Schrodinger operator H{u) for to G Vt := 
[Q-iQ+] E y i a their associated forms, 

D{f) u ) = W^T), 

f)u(f, 9) = X] K/e I fe)i2(0,l) + " fe \ ^e)x 2 (0,l)] 

ees(r) 

These self-adjoint operators correspond to the differential expression —f" + uj e -f e 
on the edges, together with the free (often called Kirchhoff) boundary conditions 
at the inner vertices, i.e. 

E /e(o)- £ f' e (o) = o (VveVnA), 

i(e)=v r(e)=v 

The coupling constants u e carry the random structure. They are picked inde- 
pendently for different edges with a probability measure /1 on 1 with supp fi = 
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[q- , q + ] . For technical reasons we have to assume that /x is Holder continuous 
with Holder exponent a and further that \i satisfies the following assumption: 
there exists r > | such that for h small 

fi([g_,g_ + /i])</ l T . (1) 

This single site measure \i defines a probability P := &) e( z E /ion fi. 

We will also need restrictions Hj^(uj) for an ebbd open A defined via the form 

D(i)L) = w^rnA) 

f)A,u,(/>#) = S [(/eUe)i2(0,l) + (^e-/ e ke)] ; 

eeE(rnA) 

which corresponds to Neumann boundary conditions at the boundary vertices 
veVndA-d. |Ku04] . 



2 The main results and the idea of their proof 

Our family of random Schrodinger operators exhibits deterministic spectrum, i.e. 
there exists a closed subset S C K s.t. a{H{u)) = S almost surely. This is a 
standard result from the theory of random operators - see, e.g., |CL90| - and 
comes from fundamental properties of our construction, especially the ergodicity 
w.r.t. lattice translations. To locate the deterministic spectrum we can consider 
the free operator H (the one with V — 0) and use some results that relate the 
spectrum of H to the spectrum of its transition operator, the Laplacian on Z d 
- see, e.g., |Ex97l ICat97| . In this way we get a(H ) = [0, oo) and hence again by 
standard theory S = oo). 

Our first claim is that in some neighborhood of inf E = g_ the operators 
exhibit pure point spectrum with exponentially decaying eigenfunctions almost 
surely: 

2.1 Theorem (Spectral/ Anderson localization) There is an e > such 
that the spectrum of H{uS) in g_ + e ] JS pure point for a.e. wfO. Further- 
more, there exists a 7 > and for each eigenfunction u associated to an energy 
in this interval a constant C u such that 

HXAiOrHI < c u ■ exp[-7d(x,0)] (x E T), 
where A x (x) is the intersection ofT with the unit cube centered at x e Z d . 

The assertion of the preceding theorem is sometimes called Anderson local- 



ization or spectral localization (see [R JLS95| for a discussion of different concepts 



of localization). An alternative and stronger concept is dynamical localization, 
see |GdB98l IDStOl] and |GK01| for more recent developments. In the context of 
our model the following result is valid. 
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2.2 Theorem (Strong dynamical localization) Let p > 2(2r — d) where r 
refers to (Qp. Then there exists an e > such that for K C T compact, each 
interval I C [E , E Q + e] and r/ G L 0O (M) with supp rj C I we have 

E{\\\X\^(H(u))xk\\}<oc, 
which in particular means that 

E{sup \\\X\>e- itH Wp I (H(u))xK\\} < oc. 
Both results will be proved by a multiscale induction as presented in detail in 



Sto01| . As the framework introduced there is general enough to include our case 
it will be sufficient to establish the necessary model-dependent estimates that are 
to be plugged into the multiscale machinery. 

For the readers convenience we will now briefly describe the idea behind the 
multiscale induction. The basic property one proves by induction is an exponen- 
tial decay estimate for the kernel of the resolvent of H^, l ^(lo). More precisely, it 
is shown that with high probability (depending on L) the resolvent of 
shows exponential off-diagonal decay. 

Note that, outside the spectrum of a Schrodinger operator, such an exponen- 
tial decay estimate is just the content of the celebrated Combes- Thomas estimate. 
We will make clear that an analogue holds for quantum graphs as well. Actually, 
this kind of argument will give the starting point of our induction procedure, the 
initial length scale estimate. More precisely, the assumption (1) on the tail of 
the single site measure implies that energies near inf E are in the resolvent set 
of H^/ L \(u) with high probability for any given L. However, keeping an interval 
near inf E fixed and letting L tend to infinity, the interval will be filled with eigen- 
values of the box Hamiltonian. Therefore the sought property, the exponential 
decay, must be deduced by a more clever argument. One important ingredient 
is the relation between resolvents of different nested boxes, cast in the form of 
a geometric resolvent identity. This will allow to conclude exponential decay on 
a large box, knowing exponential decay on smaller sub-boxes. In this induction 
step, from length L one proceeds to L a with suitable a > 1. A very important 
a priori information is necessary, the so-called Wegner estimate. Putting these 
estimates together as in [StoOlj one arrives at the desired exponential decay esti- 
mates for larger and larger boxes. To conclude, finally, that the operators H(uj) 
exhibit pure point spectrum almost surely, we need to know that the spectrum 
is indeed determined by generalized eigenf unctions. In the next section we show 
how to obtain these steps. 

2.3 Remarks (a) Our results can easily be extended to certain other cases, for 
instance, to a "rhombic" lattice, where the present method would work after ad- 
justing constants appearing in the equivalence between the Euclidean and the 
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intrinsic metric. 

(b) The results could be also extended to potentials, which are only relatively 
bounded, for instance, one can consider suitable L p (0, 1)— functions with a posi- 
tive lower bound as "single edge" potentials, following [StoOl] and numerous other 
papers; we did not take this path and treated characteristic functions as random 
potentials here exclusively for the sake of simplicity. 

(c) In a different direction, results are available for certain random quantum 
graphs, namely for random trees with random edge lengths; see the recent work 
in |ASWn61lHFn6) . 



3 The proofs 

3.1 A Combes-Thomas estimate 

The statements of this section will show how to obtain "exponential decay of the 
local resolvent" outside the spectrum. The results go back to the celebrated paper 
|CT73| and its improvement in |BCH97| . 

3.1 Theorem (Combes-Thomas estimate) Let R > 0. There exist con- 
stants c\ = ci(g_, q + , R), c 2 = c 2 (g_, q+, R), s.t. from the assumptions 

(i) AcR d T-ebdd. box, A,BcA T-ebdd., dist(A, B) =: 5 > 1, 

(ii) (r, s) C g(H^) n (-R, R) , E 6 (r, s), rj := diat(E, (r, s) c ) > 
it follows that 

\\ XA {H% - E^xbW < ci • rT 1 • e- C2 v 7 ^)*. 
Proof: Let w : A— »R be defined as w(x) := dist(x,B). By triangle inequality 

\w(y) — w(x)\ < \x — y\, 

so that || Vtolloo < 1, and this in turn implies ||«/||oo < 1 for the restriction to the 
graph. Furthermore, the functions ip(x) = e~ w ^ and (p(x) = e w ^ are uniformly 
Lipschitz continuous on all edges because 

\e w ^-e w{x) \ < sup \(expow)'(£)\-\y-x\ 
gernA 

< sup |exp(iu(0)||u/(f)l -\V-A- 
CernA 

Hence for each u G -D(h) also the functions ipu, tpu belong to D(fy), which means 
that 

\]p{u,v) := \){e- pw u,e? w v) 
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is well defined for all u, v G -D(f)). By the product rule we have the relation 

t)p(u,v) = {e-^v! | e*V) - p{{e-P w uw' \ e^v') 

-/9 2 ((e-^W | e^vw') + (3{{e- pw v! \ e pw vw') + (Vu \ v) 
= \)(u, v)-P [iuw' I v') - (u I vw')] -/3 2 (w' 2 u | v). 

" v ' 

(*) 

Referring to the term (*) above we define the symmetric form 

t(u, v) := i[{uw' | v') — (u | vw')]. 

Using 1 > m := w' 2 > one can write 

t)p(u,v) = f)(u, v) + iflt(u, v), where 
t)(u, v) = t)(u, v) — [3 2 {mu | v). 

Next we are going to show that rj^ is sectorial. From ||w'||oo < 1 one gets 

l{u) < 2\\u\\ \\u'\\ < \\u'\\ 2 + ||w|| 2 . (2) 

On the other hand, consider the operator H associated with f) and C = C(R), 
C>j3 2 + l,C>l — rm for which we have 

\\(H + C)^u\\ 2 > \\(H + P 2 + I)^u\\ 2 

= \\u'\\ 2 + ([V + (3 2 [l - m) +l]u\u) 

>o 

> iKir+n«n 2 - ( 3 ) 

It follows from (ED and (3} that 

\i(u)\<\\(H + C) l iu\\ 2 = ti + C)(u), (4) 

hence typ = t) + i/3t is sectorial and there exists an associated sectorial operator 
Hp - see, e.g., |Kat o76j. 

In the next step we are going to show the existence of a bounded operator S 

on L 2 (rn A), hsh < 1, s.t. 

t(u, v) = (S(H + C)^u\(H + C)*v) (Vu, v G D(t))). 

Let thus D{\)) be equipped with the scalar product (h + C)(-,-). By the Riesz 
representation theorem there exists a bounded operator K on D(t)) with 

t(u,v) = (t) + C)(Ku,v) 
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and by (0} we have \\K\\ < 1. Put 

S := (H + C)^K(H + c)~* : L 2 (r n A) -> L 2 (r n A). 

As + : D(fj)->L 2 (rnA) and (if + C)-* : L 2 (r nA)-> D(ty) are unitary, 
we have \\S\\ = \\K\\ < 1, and for u, v G -D(f)) we get the desired relation 

({H + C)^u\{H + C)h) = ((H + C)^Ku \ (H + C)^v) 

= t(u,v). 

Now we have to investigate invertibility of Hp — E for E G (r, s) in dependence 
on (3. Here we can use the proof of [StoOlJ (which in turn uses Lemma 3.1. from 
|BCH97| ) word by word, so we present just the result: let 



A -min^^l^l^-r) 

then for \j3\ < fa the operator T + i(3S is invertible with 

||(T + ^)- 1 ||<4^±^. (5) 

Next we will find a connection between T + if3S and Hp — E which shows that 
for |/?| < /3i the operator Hp — E is invertible too, namely 

(Hp-Ey 1 = (H + C)-*(T + i(3S)-\H + C)~i (6) 

Let / G L 2 (r n A), then 

it := (# + C)-5(T + if3S)-\H + C)-5/ g D(h) 

holds, since (i/ + C)~2 maps L 2 (A) to -D(f)). Using the definitions of T, S and it 
we can calculate for v G -D(h) the expression 

{\)p-E){u,v) = (j)-E)(u,v) + ipt(u,v) 

= (T(H + C)^u\(H + C)^v)+i(3(S(H + C)^u\(H + C)^v) 
= ((T + if3S)(H + C) l iu\(H + C)h) 
= ((H + C)- l if\(H + C)h) 
= (f\v). 

Consequently, we have u G D(Hp — E) and {Hp — E)u = /, so (jEJ) follows, and 
by (jS) we get 

\\(Hp-E)-i\\<4^. (7) 



A straightforward calculation now shows that 

{Hp - E)- l f = e Pw (H - EY l e~ l3w f, 

and therefore 

\\Xa(H - E^XbW < IIXAe-^Hoo • \\{H P - E^W • ||e^Xfl||oo- (8) 

Putting (3 := \(3\ we analyze the factors on the right-hand side. By w\b = one 
has He^xslloo < 1- The second factor is controlled by (JTJ) , and furthermore, by 
definition of f3\ there is a constant c 2 = c 2 (R) s.t. 

(3 > c 2 (R) ■ y/v(s-r). 

By assumption, w(x) = dist(x, B) > 5 for all i6i, i.e. 

IIXAe-^IU < e-f 3 - 6 < exp(-c 2 (i?) • y/rj(s - r) ■ 5). 

Combining this argument with (jSJ) we get finally the result, 

\\ X a(H - E)- X XB \\ < Cl (R) ■ tT X ■ exp(-c 2 (R) ■ yftfc - r) ■ 6). 



3.2 The initial length scale estimate 

The initial length scale estimate tells us something about the probability that an 
eigenvalue of the box hamiltonian is found inside a suitable interval. Specifically, 
we take an interval centered at the lower bound q- of the deterministic spectrum 
and we suppose that its length depends on the size / of the box. The estimate 
we are interested in will only hold for lengths larger than some initial value I*. 

3.2 Theorem (Initial length scale estimate) For each £ G (0, 2t — d) there 
exist (3 = 0{t, G (0, 2) and I* = P(r, f) such that 

F{dist{a{H^{u)),q^) < f~ 2 } < T« (9) 

holds for all T-ebdd. boxes A = Aj(0) with l>l*. 
Proof: Let 

&i,h := {cu Etl \ q e (uj) >q- + h for all e G E(T H A}. 
By the min-max principle we infer that for oj G Qi t h 

E (H£) > E Q ((—A + q^ + h)%) = q^ + h, 
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where E is the lowest eigenvalue of the respective operator. Using assumption 
(DQ) the probability of fl^h can be estimated by 



P(«j,h) > i-H^(rnA).^([g-,g- + /i]) 

> 1 - d- |A| ■ h T . 

Let £ G (0, 2r — d). Then it is always possible to choose (3 G (0, 2) such that 

f <r(2-/?)-d, 
and inserting h := Z^ -2 we get for I large 

p(nyO = i-d|A|r^- 2 ) 

= i_ rf /«-(2-/3)+^ r e 

" v ' 

<1 for i large 

> i-r*. 



3.3 The geometric resolvent inequality 

As we mentioned above, in the multiscale induction step one has to deal with 
restrictions of a Schrodinger operator to nested cubes on different length scales. 
Consequently, we need a tool that relates the resolvents of such restrictions. The 
first step on this way is the following lemma, called geometric resolvent equality. 

3.3 Lemma (Geometric resolvent equality) Let A C A' C K d be some open 
T—ebdd. boxes, H A and H A i the respective realizations of our model operator 
with Neumann b.c. Let if) G {/|rnA I / £ Cc(A)} be real-valued. Then we have 
for each z G g(H A ) n g(H A r) the relation 

R A ip = ip R A , + R A [if/ ■ D + D if/] R A ,, 

where we have denoted R A := (H A — -z) \ -^A' := (-^A' ~~ z )~ l > D is the hrst 
derivative, and all the terms are interpreted as operators on L 2 {T fl A'). 

Proof: We regard L 2 (T fl A) as a subspace of L 2 (T fl A'). In terms of the 
associated forms the assertion then reads as follows: 

(f} A -z)(i/iR AI + R A [i/-D + Di//]R At )g,w) = (if>g\w) 

(VgeL 2 (TnA),weD(t ) )); 

notice that in this case the first argument at the left-hand side, which we denote 
as u, belongs to D(H A ) and (H — z)u — if) ■ g. 

In the first step we have to show that u G ZD(fj) holds. By the product 
rule, if) e {R A >g) e G W^(0, 1) for all e G E(T H A). The continuity of if)R A >g at 
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the inner vertices of A' is clear, so the first term is controlled. Further we find 
i)'DR A , : L 2 (r n A) -> L 2 (r n A), i.e. 

R A if/DR A ,g G D(fj). 

For the analysis of the third term one has 

i>'R A , g eL 2 (rnA). 

Now R A D : W 2 (r nA)^ D(t) A ) extends to a bounded operator from L 2 (T n A) 
to D(t) A ). Indeed, we can always choose z small enough, in which case 

R((H A - z)- 1 *) = D((H A - z) 1 *) = £>(f> A ) C W 2 \T n A). 

For v G -D(I)a) we have 

IMlL(rnA) = £ 

eeE(rnA) ^° 

< ll^l|r>(f)A)> 



v1(x)dx 



i.e. D(i/A - is bounded on L 2 (r n A). Thus for ip G W 2 \T n A) and 
/ G L 2 (r n A) we get 

\{{H A -z)^'\f)\ = \(<p\D(H A -z)^f)\ 

< c ■ ||v?||i 2 (rnA) • H/IUaCrnA), 

and from here finally the boundedness of the map 

R A D = (H A - z)-*(H A - z)~*D : L 2 (r n A) — > L 2 (r n A) — > D(t) A ). 

The next step is to control the behavior of some functions at the inner vertices. 
For a fixed inner vertex of TflA let and e^^ut be the in- and outcoming edges, 
respectively, parallel to the k— th coordinate axis, and let d^{v) be the k— th 
partial derivative of the (A)— continuation of ■0. Then 

(Dip'R A >g | w)l 2 (iTiA) = £ (Dip' e R A ,g e \ w e ) L2(0 ,i) 

eee(rnA) 

ee£(rnA) 
= -(i/j'Ra'9 I w'k 2 (rnA) 

+ £ Ed fe ^){(i?A'<He fc , in (l) " (i?A^^)e fc . out (0)} 

■u inn. vertex fc=l v 

=0 by continuity at inner vertices 

= -(tfj'R A >g\ w') L2 ( TnA) . (10) 



10 



The following calculation now finishes the proof: 

(f) A -z)(u,w) = (t) A -z)(tjjR A ,g,w) + ((ip'-DR A ,+DiP'R A ,)g\w) 

® (ty R A ,g)' | w 1 ) + ((V - z)1> R A >g \ w) 
+ (i/j'(R Ag y\w)-(^R A g\w') 
W R A <g | w') + ^(R A g)' \ w') + (ip'(R A g)' \ w) 
+((V - z)if> R A >g | w) - (il>'R A g \ w 1 ) 

((R A , g y\(i>w)') + ((V-z) R A ,g\i>w) 
(h A ' - z)(R A g,tpw) 
(g | ipw) 
(ipg\w). 



tfi real val. 



We will next prove another preparatory lemma after which we will be ready 
to state the main theorem of this section. 

3.4 Lemma Let Q C Q, C M d be a T-ebdd. domains, dist(dtt, dQ) > 0, E G R 
and g e L 2 (TnQ). Then there exists C = C(g_, q+, E) s.t. for all u G W^TnQ) 
with 

W I y')i 2 (rnn) + (Vtt | <^)L 2 (rnn) = I ^(mo) (Vp G n fi)) 

it izoids that 

IKIIz 2 (mn) < C(lk|U 2 (rnn) + IMk(mn))- 
Proof: By construction, dist(5f2, dCt) > 1, hence there exists a vector 

if> G {f\mn I / G C c (fi), supp /' C {x G | dist(x, Q) < 1}} 

with < if) < 1, V'lrnn = 1 an d ll^'lloo < C(d). Let w := m/; 2 , then w G 
VF 2 X o(r n Q), and by product rule we find 

(u' | w')i 2 (rnn) = (ipu' | ^w') + 2(ipu' \ uip'). 

Using V := V — E and support properties of the functions involved we get 

\\iju'\\ 2 = («' I w') - 2{^u' | ui//) 

= (g | w) — (Vu | w) — 2(ipu' | 

< NIN| + |(^I^)I + 2|MIIHI||^IU 

< \\g\\\\u\\ + C( q ^q + ,E)\\u\\ 2 + 2C\\i;u'\\\\u\\. 
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We consider the latter as a quadratic inequality in || ||, and find after some 
simple manipulations, that it can only be fulfilled for 



< yc 2 + c\\u\\+ . \\ g \\ 

2VC 2 + C 
= C(q„,q + ,E)(\\u\\ + \\g\\) 

By "0lrnc = 1 the assertion follows. □ 

Before we come to the main point we introduce some notation. A T— ebdd. 
box A = A L (x) is called suitable, if x G Z d , L G 6N \ 12N and L > 42. For such 
boxes we define 

A int (z) = A LM{x) := A L/3 (x), 
A out (a;) = A i)0ut(a: ) := A L (x) \ A L _ 12 (x) 

and write for the respective characteristic functions on the graph: 

v int int int -v out -\/ out •v out 

AA — XA L (x) -~ ArnA int (x)' aa — AA L (x) -~ ArnA int (x)- 

In general the symbol xa for a T— ebdd. domain is to be understood as XvnA- 

3.5 Theorem (Geometric resolvent inequality) Let A C A' C M. d be suit- 
able T—ebdd. boxes. Let further A C A mt and B C A' \ A be T—ebdd. domains, 
/od bounded and E G I . Then there exists C geom = C geom (g_, q + , E) s.t. 

\\XbRa>{E)xa\\ < C geom -\\xBRA'(E)xT\\ \\xTMe)xa\\- 

Proof: Let x G Z d be the center of A. We choose ip G {/|rnA | / G C C °°(A)} 
real-valued with supp / C A L _ 4 (x) s.t. <p = 1 on A L _ 8 (x). This can be certainly 
achieved, with ||y?'||oo bounded independent on A. 

Let fl := intA out , i.e. dist (dfl, supptp') > 2. By the geometric resolvent 
equality (Lemma 13.31 ) we have 

IIxb-Ra'XaII = ||xa-Ra'Xb|| 

= \\xa(<pRa> - Ra<p)xb\\ (<p\a = 1,<p\b = 0) 

Lem = ° \\xa(vRaW + ¥ >>D)R a ,xb\\ 

< WxavRaDip'Ra'XbW + \\xa<pRa<p'DRa'Xb\\ ■ 



(*) (**) 

We start with the analysis of (*). If fl := int (Ai_ 2 (x) \ Al- W (x)) it holds that 

(*) = \\XA<pRADXnXnLp'RA'XB\\ 

< \\<p'\\oo \\xa^RaDxq\\ WxsiRa'XbW- 

y v ' 

(***) 

12 



The term (***) can be now controlled with the help of Lemma [3T4l We put 

/ G L 2 (T fl A), g:=XAf, u:=R A g. 
Then u G £>(h) and 

(fj A - E)(u,w) = (g\w) 

for all u» G D(1)a). Furthermore, we have g\n = as well as dist (dfl,dfl) = 1. 
Consequently, Lemma I3T41 is applicable and it gives 

||Xn u 'll < C l i(g_,g+,J) ||u||i 2 (rnn) 

= Cx(q_,q + ,I) WxnRAXAfW, 

i.e. 

(***)< d g+, J) ||xa ^aXaII- 
The term (**) can be treated in a similar way. □ 

3.4 The Wegner estimate 

The Wegner estimate represents a statement about the probability that the op- 
erator Hj^(u>), restricted to a T— ebdd. box A = Ai(x) centered at x G Z d , will 
have an eigenvalue near some fixed energy. Typically - and sufficiently for our 
multiscale analysis - this probability is polynomially bounded in terms of the box 
volume. 

3.6 Theorem (Wegner estimate) For each R > there exists a constant Cr 
such that for all T—ebdd. boxes A = Aj(i), i G Z d , and all intervals I C (-R, R) 
of length |/| the following estimate holds: 

F{a(H%(u))Dl^<J)} < C R -\A\ 2 -\I\ a . 

Before we start with the proof let us recall the following elementary lemma from 
[StoOOj . 

3.7 Lemma Let J be a hnite index set, \x a Holder continuous probability mea- 
sure on M. d with Holder exponent a, /i J := ®j G jAi the product measure on M. J . 
Let $ : M J — > R a monotone function, for which there are constants 5 and a > 
s.t. for all t G [0, 5], q G R J we have 

$(g + *(l,...,l))-$(g) >t-a. (11) 

Then for each interval I of length smaller than e < aS the following estimate 
holds: Q 

//({g:$(g)G/})<|J|-(£) . 
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Proof of Theorem \3.(k We start with an estimate for the number of eigenvalues 
smaller than a given energy R. To this aim we define the Neumann-decoupled 
operator _A^' dec via its quadratic form 

t)f de Vi9) ■= E (/VWnA). 

ee-E(rnA) 

By a direct calculation the eigenvalues of this operator are ^n 2 , n G No, with 
the multiplicity §{E(T fl A)} < d ■ l d = d\A\. Hence there exists a constant Cr 
s.t. for the n-th eigenvalue, counting multiplicity, it holds that 

E n {-Af dec ) > R forn > C R \A\. 

Now we have 

H»(co)>(H + q ^>-Af d - 

since g_ > by assumption, and thus by min-max principle the corresponding 
inequality for the n-th eigenvalues. Using the previous inequality we get 

PH^H)n/^0}< £ F{E n (H»(u;)) G I}. (12) 

n<C R -\A\ 

Next we estimate the terms of the sum by means of Lemma [3771 Because of the 
independence of H a (uj) of coupling constants outside A we have 

F{E n (H A (u)) G /} = fi E{r) {u\E n (H A (u;))eI} 

= ti E{rnA) & = (w e ) ee£( rnA) I E n (H A (u)) G /}. 

By := E n (H A (u)) = E n (H A (u)) a monotone function on R E (T n A) is 

defined, and it fulfills condition ([Til because 

H A (u + t(l,...,l)) = -A+ J2 (^ + *)Xe 

esE(rnA) 
= H A (u)+t. 

Hence by Lemma 13771 we have 

F{E n (H A (uj)) g /} < p(rnA)-|/r 

< d\A\\I\ a , 

which in combination with ( fT2l ) yields the assertion. □ 
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3.5 Expansion in generalized eigenfunctions 

Now we come to the last statement needed for the multiscale analysis, namely 
that polynomially bounded generalized eigenfunctions exist spectrally a.s. 

We want to use the main result from [BMSt03| . that gives the polynomial 
boundedness in terms of the intrinsic metric (see |Stu94| ) generated by the free 
Laplacian H on the graph. Using the embedding of our graph into M. d it can 
easily be seen that the intrinsic metric is equivalent to the Euclidean one on IR d , 
and consequently, after adjusting some constants the statement can be written 
in terms of absolute values as well. We start by checking the assumptions of 
|BMSt03| . First of all one has to show that the form f)o associated with the free 
Laplacian is a Dirichlet form. Note that || • ||^ is equivalent to the norm || • Hw^r) 
so f)o is closed. For u G -D(f)o) which is real- valued we have |m| G D(t) ), and 
therefore ^ 

f)o(M) = ^2 / ( s S nu e(x)u e (x)) 2 dx = f)o(w)- 
eeE(r) ^° 

If u is in addition nonnegative, we have u A 1 G -D(f)o) and 

fjo(«Al)= ^2 / u 'e( x ) 2 ■ l[ Ue <i](x)dx < ffo(lt). 



Obviously f) is strongly local and regular - see, e.g., [BMSt03j for definitions. 

The next point is that the volume of balls with respect to the intrinsic metric 
q does not grow too fast as R^-oo. Because the graph is embedded into M. d and 
the intrinsic metric and the || • ||i-metric are equivalent, the volume of the ball 
B e R (x) can be estimated by the number of edges contained inside a box A 2 r(x) 
and hence by CdR d for large R. 

Finally, the third assumption to be checked is that e~ tH ° is bounded as a 
map from L 2 {T) to Loo(r) for some t > 0. To this aim we employ the following 
extension of the ultracontractivity result [KMS06], Lemma 3.2, demonstrated by 



using the same method as in the cited paper. 
3.8 Lemma For t G (0, 1] it holds that 

||e-^|| L2(rHW r) < cH. 

Proof: By |Ou05| . Thm. 6.3 ff, see also [Na58l lFSt86l IDav89j . it is sufficient to 
show that 

U 2 (r)<C||/||M|/||I i(r) 



for / G D{\)) n Li(r). Now, by [Ga59l INi59] . or by |Ma85l Sect. 1.4.8], we have 
the following Nash type inequality for u G W^i^O, 1): 



MU 2 (o,i) < ci • (||w'|U 2 (o,i) + IMUi(o,i)) 3 • ||«| 
< c i ' IpIIwj^o.i) ' IMIii(o,i)' 



L!(0,1) 
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where in the second step the Holder inequality has been applied to u • 1. For 
/ G D{bt) n Li(T) we have by another application of Holder inequality 

\\J\\L 2 (T) ~ E « 

L 2 (0,1) 

eG-E(r) 

- C l E ll/e||^ri( ,l) ' H/ e llLi(0,l) 
ee£(r) 

< 4( E H/.ll^(0.1)) 4 • ( E H/ellw)) 1 

eeE(r) eeB(r) 




With these assumptions, given using the arguments in the opening of the 
section, |BMSt03) yields the following result: 

3.9 Theorem For spectrally a.a. E e cr(H) there exists a generalized eigen- 
function tp 

(i + |.| 2 )-T^ Gj L 2 (r). 

satisfying for any m > ^±1 , 

This completes the necessary input for the use of Theorem 3.2.2 from [Sto()f| and 
thus the proof of Theorems 2.1 and 2.2. 
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